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In this note we address the problem of solving for the positronium mass spectrum. We 
use front-form dynamics together with the method of flow equations. For a special choice 
of the similarity function, the calculations can be simplified by analytically integrating 
over the azimuthal angle. One obtains an effective Hamiltonian and we solve numerically 
for its spectrum. Comparing our results with different approaches we find encouraging 
properties concerning the cutoff dependence of the results. 

1. Introduction 

Solving for QCD bound states from first principles remains a yet unsolved 
problem. There is some hope of simplifications in the framework of front-form 
dynamicsEI, because of a simpler vacuum state. Different methods have been devel¬ 
oped over the last years to cope with the stunning problem of constructing an effec¬ 
tive Hamiltonian out of the ’infinite’ canonical QCD Hamiltonian^ The method of 
flow equationsB is closely related to the similarity transformations 0 . Both methods 
are based on unitary transformations to block- or band-diagonalize the Hamiltonian. 
So far no calculations using the fully relativistic andxovariant effective light-cone 
Hamiltonian were performed in the manner of Refs, uo, where first the many body 
part of .the problem is resolvedu. We present numerical results in the similarity flow 
schemecl applied to positroniuim with a special choice of the similarity function. We 
compare to the results of Ref. □, and point out the specific virtues of the methods. 

2. Flow Equations and Similarity Renormalization 

Both the flow equations of Wegner @ and the similarity renormalization of Wil¬ 
son and GlazekS can be retrieved as special cases of the more general similarity 
flow framework E9. The argument goes as follows. We construct an effective Hamil¬ 
tonian H' from the canonical Hamiltonian H (regularized at some scale A) by a 
similarity transformation, H' = UHU'. The generator of the transformation is 
anti-hermitian, rf = —r). The Hamiltonian is considered a function of the flow pa¬ 
rameter l, with the bare and effective Hamiltonians being H(l = 0) and H(l —> 00), 
respectively. Its change with respect to the flow parameter is given by the commu¬ 
tator with the generator 77 of these unitary transformations 

^r 1 = m,H(i)}. (l) 


1 
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The goal is to choose an 77 such that the transformed Hamiltonian has the form 

FjW- ( 2 ) 

The flow parameter is connected to the scale A by l = 1 /A 2 , and the free particle 
energies Ei are defined by the free Hamiltonian 

H d 0) \i) = Ei\i). ( 3 ) 

If the so-called similarity function /(A) vanishes for vanishing A, then will 

become a (block) diagonal operator in this limit. The generator 77 of the uni¬ 
tary transformations is by construction the commutator of the diagonal part of the 
Hamiltonian with its complement, 77 = [Hd,V\. The two schemes of Wegner and 
Wilson-Glazek are defined by different choices of the generator 77. Note that the 
similarity function /(A) is still arbitrary in both schemes. 

3. Positronium Model 

A solution of Eq. (JI) was described in detail in i, so we can be brief here. 
One expands the Hamiltonian and the generator 77 into a power series in the bare 
coupling g. With the energies Ei depending on the flow parameter only in second 
order, we can solve the differential equation for the Hamiltonian, Eq. [I] in the energy 
basis defined by Eq. (j|) order by order. Using the definition 

V%\l)=f{l)vS\l), ( 4 ) 

suggested by Eq. we get the solutions in the second order 

Hffil) = Hg\l = 0)+ ( 5 ) 

Jo 

= \f\l = 0)+ f ( 6 ) 

Jo 

where the superscripts ( d ) and ( V ) denote the diagonal and the particle number 
changing part, respectively. We are supposed to take the bare cutoff A, which 
defines the terms at l = 0, to infinity. 

To evaluate Eqs. (||) and (](]) and to obtain the matrix elements of the Hamilto¬ 
nian, we choose the particle number conserving part of the Hamiltonian as diagonal, 
and we reduce the particle number violating blocks of the Hamiltonian to zero with 
the flow equations. This solves en passant also the multi-particle problem. Instead 
of having to truncate the Fock space a la Tamm-Dancoff, we are now dealing with 
isolated blocks of definite particle number. The effective matrix elements are thus 
obtained by the process 


Ueff = lim (V 9en + V PT ) . 

A—>0 v ' 


( 7 ) 
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The interaction can be read off the structure of Eqs. (|) and (§): the part without 
the integral ( V PT ) is the one obtain by usual perturbation theory. The second one 
(V gen ) is generated by the flow of the Hamiltonian. In other words, by reducing 
the off-diagonal matrix elements, we are inducing changes on the diagonal. 

To calculate the matrix elements, one evaluates the associated diagrams, ap¬ 
plying light-cone perturbation theory liiJ. The electron [positron] momenta and 

currents are l£ = ( k' e — k e Y[lg = {kg — fcg) M ], and j(l e ) M = u{k' e ) r y IJ 'u{k e )[j{lg) IJ ’ = 

/ 

u(k' g ) r y^ u (k g )\, respectively. Here, kg and kg' are the electron momenta before and 
after the interaction. For the numerical calculations we use relative coordinates: 
x = -pq- is the longitudinal momentum fraction, and k± is the transverse momen¬ 
tum. We obtain 


ygen 


V\ 


PT 


J (leTMle) f 
Ai J X 


dfy{ AQ 
d\' 


f\'{^2)dX' + (1 2 ) 


V 


fy{A 1 )fy{A 2 ) + 3+ ^ 3+ ^ (T* - u). (8) 


Here, T* = | (k e + kg) 2 + \{k' e + k' s ) 2 is the average kinetic energy before and after 
the interaction and w is one of the (unknown) eigenvalues of the full Hamiltonian. 
The latter ambiguity is no problem in the formalism considered here, because the 
perturbative term vanishes when the scale A goes to zero. The energy denominator 
is given by 

V=\x-x'\(T* ( 9 ) 

and we used the definition A j = Ai/(x — x'), where 


Aq = mj 


(x — x ') 2 
xx' 



2 kj_k' ± cos (<p — ip'), 


( 10 ) 


and A'2 = Ai(x —> 1 — x, x' —> l—x r ). The similarity function f\(A) in the inter¬ 
action is still at our disposal. In the electron-positron sector we have the integral 
equation 


0 = 


+ 


^ J dx l d 2 k , j_{x,k ± -Xi,\ 2 \V e s\x',k , j_-, 

Aq ,X 0 


(ii) 


^2 )’ t l’n{x', k'j_; Aq, A2). 


The cutoff enters into the problem via the definition of the integration domain. It is 
restricted by the Lepage-Brodsky cutoff on the kinetic energies (m 2 -|-fc^)/a:(l — x) < 
A 2 + 4 m 2 . To simplify the numerical calculations we integrate out the azimuthal 
angle in the problem (&_l = kj_e' ltp ) by substituting it by the discrete quantum num¬ 
ber J z , cf. Refs. Bta. To be able to perform the corresponding integral analytically, 
we have to choose a special similarity function 


/a (A) = 9(X 2 — |A|). 


( 12 ) 
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Figure 1: Positronium spectrum(a = 0.3, A = 1.0 m, N\ = N 2 = 21): bound state 
masses in units of the electron mass in different J z sectors. 

This sharp cutoff leads to difficulties in the calculations, namely the collinear singu¬ 
larity is not canceled exactly anymore. The scale integrals, however, become very 
simple in the limit A —■> 0, yielding a theta function, 0(A, — A j). 

The actual matrix elements, integrated analytically over the azimuthal angle 
differ only by a re-definition of the spin-dependent function Int(n) from the ma¬ 
trix elements of the Pauli ansatz, cf. e.g. Ref. 0 . This makes it possible to use the 
numerical techniques and the computer code of Ref. 0 to calculate the positronium 
spectrum. The matrix elements used in the present note have an additional sin¬ 
gular spin-independent part. We argue that this is an artifact of the choice of the 
similarity function and will omit it in the numerical calculations. We comment on 
the justification of this step below. The interaction used here is the same as listed 
in Appx. F of Ref. 0 , with the following changes. 


Gi(x,k±;x\k' ± ) - 

-*■ G\ + — T^Intdl - n|)< 5 nj i 

[X — x'Y 

G 2 (x,k±-,x',k' x ) - 

-* G 2 +\ jlnt(|n|)<5 nj o, 

(x — x’Y 


and in all Gi the function Int(n) is replaced by 

Int(n) = 9(ai — a 2 )Int\{n) + 0{a 2 — ai)Int 2 {n ), (13) 

where 

Oj(a? — 4 fc^fcj 2 ) -1 / 2 — 1 
2k±k l J _ 



Inti(n) : 


(14) 
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Figure 2: Cutoff dependence: (a) Eigenvalues. Below: triplet (upper curve) and 
singlet (lower curve) ground state. Above: first excited states (n=2). (b) Hyper- 
fine structure coefficient C^f- The cutoff is given in units of the electron mass 
(a=0.3,JVi=25,JV2=21). 




Here 


ai 


(x — x') 2 
xx' 


+ k 2 ± + k'l 




(15) 


and(i 2 = o,\{x —> l—x,x' —>• 1 — x'). Note ilrat the analogue of the latter expressions 
in the calculations with the Pauli ansatz □, is their average a = + a2 )/2. 


4. Numerical results 

The Hamiltonian matrix elements were derived by applying the flow equation 
scheme to the positronium problem in front-form dynamics. Contrary to preceding 
work, at this point of the calculations we use a non-perturbative method to extract 
the spectrum of this Hamiltonian rather than light-cone bound state perturbation 
theory. We have to solve the eigenvalue problem Hlc\ti) = M 2 \n) or equivalently, 
the integral equation, Eq. ©■ We use the algorithm set up in Ref. 13 . For details 
of the calculations and numerical methods applied, see there. To be able to trace 
possible violations of rotational invariance (a non-trivial issue in front-form dynam¬ 
ics), we chose to work with an unphysically large coupling constant a = 0.3. This 
is possible, because the integral equation is an algebraic function of the coupling. 

The results of the computations are compiled in Fig. [lj. We get the expected 
Bohr spectrum, the hyperfine splitting and even the correct multiplet structure. 
The multiplets are (almost) degenerate. We find exponential convergence of the 
eigenvalues when approaching the continuum limit (N —> c«), as in Ref. u. 

The eigenvalues depend weakly on the cutoff A, cf. Fig. ||(a), except for the 
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n 

M 2 n [J z = 0) 

M 2 pt 

AM 2 

n 

M 2 (J z =0) 

M 2 pt 

AM 2 

i 

3.906908 

3.900002 

0.6906 

10 

3.990189 

3.989675 

0.0513 

2 

3.915533 

3.910673 

0.4860 

11 

3.990239 

3.989825 

0.0413 

3 

3.976824 

3.975860 

0.0964 

12 

3.990345 

3.989875 

0.0470 

4 

3.977998 

3.976533 

0.1465 

13 

3.990368 

3.989875 

0.0492 

5 

3.978166 

3.977037 

0.1129 

14 

3.990423 

3.989945 

0.0478 

6 

3.978694 

3.977206 

0.1488 

15 

3.990432 

3.989945 

0.0487 

7 

3.978705 

3.977206 

0.1500 

16 

3.990441 

3.989955 

0.0486 

8 

3.978936 

3.977441 

0.1495 

17 

3.990506 

3.989955 

0.0551 

9 

3.989893 

3.989476 

0.0417 

18 

3.990519 

3.989984 

0.0536 


Table 1: Positronium spectrum (a=0.3, A=50 m, Ni=N 2 =21): present results 
{M 2 ), perturbation theory ( M PT ), and difference (A M 2 =M 2 —M PT ) in percent. 


singlets. From analytic arguments, we expect the eigenvalues to diverge logarithmi¬ 
cally because we used flow equations derived up second order only. However, as we 
shall see, the coefficient of the hyperfine splitting is very well described at moderate 
values of the cutoff. We therefore compare our eigenvalues to results of equal time 
perturbation theory at A = 50m/. The eigenvalues are listed in Table [3] and agree 
cum grano salis with the known results. Note, however, that perturbation theory 
might be not very reliable at this large coupling. 

A test for the theory is the value of the coefficient of the hyperfine splitting, 
C h f = ( Mtripiet - M singlet) /ma^. The plot of C h f versus the cutoff, Fig. |(b), 
is encouraging: we obtain a smooth curve, converging for large cutoffs to a value 
Chf 0.2825. From equal time perturbation theory we would expect a value 

between 1 /3 (order a 4 ) and 0.2379 (order a 6 log a), which is exactly what we obtain. 

Let us focus on violations of rotational invariance, cf. Ref. 0. We fitted the 
cutoff dependence of the discrepancy of corresponding eigenvalues of different J z to 
a polynomial in logA/m/. The (negative) linear coefficient is smaller than 10 —3 , 
and the constant term is ca. 10 -4 . This means that the discrepancy between the 
triplet levels is 1% of the relevant (hyperfine splitting) scale at a cutoff of one 
fermion mass, and rises to roughly 10% when A = 18m/. This suggests that to 
obtain full rotational invariance, one has to go to higher orders in the derivation of 
the flow equations. 

The original idea of the flow-equation approach implies the use of a smooth 
similarity function to cancel the collinear singularity completely. We had to omit 
the spin independent singular part by hand to be able to perform the integration 
of the numerical counterterms. We argue that the singular part is an artifact of 
the choice of the similarity function. Indeed, it can be shown that the singular part 
of the matrix elements vanishes for the standard flow equation similarity function. 
There is a danger that this choice might change the eigenfunctions significantly. We 
found, however, that the wavefunctions are almost identical with those of Ref. □. 
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5. Discussion 

We presented the full positronium spectrum and wavefunctions using flow equa¬ 
tion techniques in all sectors of the angular momentum J z . The results are encour¬ 
aging especially when looking at their cutoff dependence. Quite in general we find 
a weak, logarithmic, smooth dependence on the cutoff. The comparison to results 
of equal time perturbation theory gives good agreement. With respect to its cutoff 
related properties, the flow equation scheme seems to work better than the model of 
Ref. u, inspired by the method of iterated resolvents. This is expected by construc¬ 
tion of the methods and new ideas have been proposed to deal with renormalization 
issues in the latter method i. The inclusion of the annihilation channel would be 
a straightforward calculation, cf. E3. We found that rotational invariance is obeyed 
to a large extent on the numerical level and withstood therefore from implementing 
this channel here. It remains to be investigated how the singular terms influence 
the spectrum. We argued that they are relics of the chosen similarity function and 
omitted them. The numerical results seem to support this claim. In general, it 
would be interesting to calculate the effective Hamiltonian to a higher order in the 
bare coupling constant. This would be useful to find out about the structure of the 
generated (irrelevant) operators and also to test the improved cutoff dependence of 
the results. 

Acknowledgments: This work was done while the author was a Minerva fellow 
at the Weizmann Institute, Israel. The author is indebted to E. Gubankova for 
enlightening discussions about her analytic calculations. 
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